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摘要 : 不 同 的 品种 抗 性 不 同 ， 为 进一步 探究 不 同 火龙 果品 种 之 间 的 抗 性 差异 ， 为 后 续 火 龙 果 
抗 性 育种 提供 参考 ， 该 研究 利用 lumina HiSeq 2000 测序 平台 对 ‘普通 白肉 ' BR) M SJEA 
ZRA EY) 两 个 品种 进行 转录 组 测序 分 析 ， 并 参考 GO Ontology, KEGG 等 公共 数据 
库 对 差异 表达 基因 进行 功能 分 类 与 富 集 分 析 。 结 果 显 示 : (1) BR 与 EY 共有 14248 个 差 
异 基因 ， 其 中 5 446 个 基因 上 调 ，8 802 个 基因 下 调 。(2) 相关 GO 功能 分 析 表 明 这 些 差异 
基因 主要 参与 酶 催化 活性 、 细 胞 组 分 、 代 谢 过 程 等 ， 其 中 参与 氧化 还 原 酶 活性 的 349 个 差异 
基因 在 BR 中 表达 量 上 调 。 (3) KEGG 通路 分 析 显 示 ， 大 部 分 差异 基因 富 集 在 新 陈 代谢 、 
生物 合成 等 , 其 中 参与 角质 、 木 栓 质 和 蜡 质 生物 合成 的 差异 基因 有 12 个 , 如 CYP86 与 CERI 
等 。 参 与 氧化 还 原 酶 活性 的 差异 基因 在 BR 中 较 EY 表达 量 上 调 ， 且 显著 富 集 ， 表 明 BR 与 
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EY 在 生长 发 育 、 细 胞 代谢 过 程 差 异 显 著 。 参 与 角质 、 木 栓 质 和 蜡 质 生物 合成 的 差异 基因 在 
BR 中 表达 量 上 调 ， 此 类 基因 在 BR 中 具有 较 高 的 表达 量 且 显著 调集 ， 表 明 BR 可 能 具有 较 
强 的 抗旱 和 抗 病 能 力 。 
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Abstract: Different varieties have different resistances. In order to further explore the resistance 
differences in different varieties of pitaya fruit and provide a reference for further study on 
breeding of pitaya fruit resistance. Our study using Illumina HiSeq 2000 sequencing platform to 
sequencing the transcriptome of ‘Putongbairou’” (BR) and ‘Ecuador Yellow? (EY). Functional 
classification and enrichment analysis of differentially expressed genes (DEGs) were performed 
by reference to GO Ontology, KEGG and others databases. The results were as follows: (1) There 


收 稿 日 期 : 2021-07-04 

基金 项 目 : 国家 重点 研发 计划 项 目 (2019YFC0507503, 2019YFC0507504-03); 广西 创新 驱动 发 展 专项 资 
& EF} AA2016100403-01); 广西 自然 科学 基金 (2019GXNSFBA245097); 广西 喀斯特 植物 保育 与 恢复 生 
态 学 重点 实验 室 基金 (19-050-6) [Supported by National key Research and development program 
(2019YFC0507503, 2019YFC0507504-03) ;  Innovation-Driven Development Special Fund of Guangxi 
(AA2016100403-012; Natural Science Foundation of Guangxi (2019GXNSFBA245097); Fund of Guangxi Key 
Laboratory of Plant Conservation and Restoration Ecology in Karst Terrains (19-050-6)]. 

作者 简介 : FEE (1992-), 硕士 , 助理 研究 员 , 主要 从 事 植 被 生态 学 研究 , (E-mail) jxlee14082082163.com. 
”通信 作者 : 陆 树 华 ， 副 研究 员 ， 主 要 从 事 恢复 生态 学 研究 ，(E-maiD lushuhual3@163.com。 


yx 


were 14 248 DEGs between BR and EY, of which 5 446 genes were up-regulated and 8 802 genes 
were down-regulated. (2) GO functional analysis showed that these DEGs are mainly involved in 
enzyme catalytic activity, cell components, metabolic processes, etc. Among them, there were 349 
differential genes involved in oxidoreductase activity. (3) KEGG pathway analysis showed that 
most of the DEGs were enriched in metabolism, biosynthesis, etc. , and 12 key genes such as 
CYP86 and CERI involving in Cutin, suberine and wax biosynthesis. We found the expression of 
DEGs involved in oxidoreductase activity were higher in BR than that in EY, which significantly 
enriched, indicating that may be differences in growth and cell metabolism between BR and EY. 
DEGs involved in the biosynthesis of Cutin, suberine and wax biosynthesis were up-regulated in 
BR. Such genes have higher expression in BR, and significantly enriched, suggesting that BR may 
be superior drought and disease resistance than EY. 
Key words: pitaya fruit, transcriptomic, different varieties, differentially expressed genes 

火龙 果 属 仙人 掌 科 (Cactaceae )、 量 天 尺 属 (Hylocereus Britton & Rose) EX rz HE Js 

(Selenicereus (Berger) BritIsraelton and Rose) 植物 ， 营养 美味 以 及 独特 的 外 观 使 其 成 为 消费 

者 追捧 的 新 兴 热 带 水 果 〔 陶 金 etal, 2014)。 近 几 年 ， 我 国 火龙 果 产 业 发 展 迅 猛 ， 据 最 新 数 
H- WUE 2020 年 底 ， 全 国 火 龙 果 种 植 面积 已 超过 100 万 亩 ， 并 形成 以 广西 、 广 东 、 海 南 、 
云南 和 贵州 五 省 (区 ) 为 核心 的 主 产 区 。 特 别 近 年 来 在 广西 石 漠 化 区 推广 种 植 成 功 ， 形 成 贫 
困 山 区 脱贫 致富 双赢 的 “ 果 化 模式 "火龙 果 耐 晶 、 耐 将， 有 利于 石山 地 区 水 土 保持 ， 在 改善 
生态 环境 的 同时 又 可 增加 可 观 的 经 济 效益 。 作 为 一 种 新 兴 水 果 , 近年 来 国内 外 对 火龙 果 的 研 
究 也 越 来 越 深入 , 主要 集中 在 不 同 品种 的 理化 性 质 和 营养 价值 、 功 能 性 物质 提取 、 和 栽培 技 术 、 
病虫害 研究 〈 赵 志平 和 杨 春霞 ，2006; Menezesetal,2015; 陈 煜 etal,2017; 924516, 
2017; Kee etal., 2019; Magalhães et al., 2019; Rahmati et aL, 2019)， 转 录 组 测序 与 微 卫星 
引物 开发 〈 杨 仕 美 等 , 2018; 武志 江 等 , 2020), 属 间 亲 缘 关 系 、 杂 交 及 新 品种 的 培育 (Tel-Zur 
et al., 2004a; Tel-Zur et al., 2004b; Tel-Zur et al., 2011; Ortiz-Hernández & Carrillo-Salazar , 
2012; Tel-Zur et al., 2012) 等 方面 。 
量 天 尺 属 包含 有 14 个 种 ， 世 界 范 围 内 广泛 栽培 的 火龙 果 主 要 为 量 天 斥 属 植物 ， 有 H. 
undatus. H. monacanthus 和 H. megalanthus COrtiz-Hernández & Carrillo-Salazar, 2012). 。 通 过 
杂交 ， 目 前 火龙 果品 种 繁多 ， 文 献 有 报导 收集 218 份 种 质 资源 〈 黄 凤 珠 等 ，2019)， 也 有 相 
关 文 献 报导 ,在 种 质 资源 保存 过 程 中 发 现 , 不同 品种 的 火龙 果 其 抗旱 和 抗 病 能 力 差异 显著 ( 梁 
秋 玲 等 ，2011; 李 润 唐 等 ，2017)。 目前 通过 RNA-seq 技术 ， 可 获得 大 量 的 转录 组 序列 信息 ， 
开展 基因 表达 与 功能 分 析 , 构建 信号 通路 等 。 本 文 对 火龙 果品 种 ‘普通 白肉 "和 “厄瓜多尔 黄龙 
植株 幼 嫩 茎 段 组 织 进行 转录 组 序列 从 头 组 装 ， 通 过 Unigene 表达 丰 度 分 析 、 差 异 表达 分 析 、 
GO 和 KEGG 富 集 分 析 等 ， 从 转录 组 层面 阐明 火龙 果 不 同 品种 之 间 的 差异 ,为 火龙 果 产 业 的 
发 展 、 推 广 应 用 提供 一 定 的 理论 支持 。 
1 材料 与 方法 
1.1 材料 

火龙 果品 种 ‘普通 白肉 ' (BR， 以 下 简称 “白肉 " ) “厄瓜多尔 黄龙 ' CEY， 以 下 简称 “ 黄 
JE') 均 取 自 广 西 壮 族 自治 区 中 国 科学 院 广西 植物 研究 所 火龙 果 研 究 基 地 【图 1)， 选 取 标 记 
这 两 个 品种 茎 段 上 开始 出 芽 的 位 置 , 待 其 生长 7d, 取 幼 嫩 茎 段 , 每 个 品种 设 3 个 生物 学 重复 ， 
采集 后 立即 放 入 液 氮 中 ， 交 由 北京 诺 不 致 源 生物 信息 科技 有 限 公 司 进行 测序 。 
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A. “普通 白肉 (BR) B. “厄瓜多尔 黄龙 (EY) 
A. *Putongbairou (BR) B. ‘Ecuador Yellow’? (EY) 


图 1 试验 用 火龙 果品 种 
Fig. 1 Pitaya fruit was used in the experiment 
1.2 文库 构建 
用 带 有 Oligo (dT) 的 磁 珠 富 集 火龙 果 样 品 中 的 mRNA。 随 后 将 打 断 的 mRNA 作为 模板 ， 

合成 二 链 cDNA， 再 进行 纯化 ， 末 端 修复 、 加 A 尾 并 连接 测序 接头 。 最 后 进行 PCR 扩 增 ， 纯 
化 PCR 产 物 ， 构 建文 库 。 待 文库 检测 合格 后 上 机 测序 ， 采 用 Inlumina HiSeq»2000 PE150 对 样 
品 转录 组 文库 进行 高 通 量 测序 。 
1.3 De novo 拼接 和 注释 

过 滤 测 序 得 到 Raw read， 获 得 的 Clean read。 在 Trinity 软 件 (Grabherr et al., 2011) 中 进 
行 De novo 拼 接 。 在 NCBI blast 中 将 拼接 完成 的 Unigene 与 NCBI 非 宛 余 和 蛋白 序列 数据 库 、 NCBI 
核酸 序列 数据 库 、Swiss-Prot 数 据 库 、KOG (Clusters of orthologous groups for 
eukaryoticcomplete genomes) 数据 库 和 KEGG (Kyotoencyclopedia of genes and genomes) 数 
据 库 进行 比 对 注释 。 利 用 HMMER 3.0 package Protein family (Pfam) 进行 比 对 Chmmscan 
e-valuex1*?) 。 在 Blast2go (http:/www.blast2go.comy/ b2ghome) 软件 和 WEGO 软 件 中 进行 GO 

(Gene Ontology) 注释 、 分 类 统计 。 

1.4 差异 基因 表达 分 析 

采用 DESeq2 (Love et al, 2014) 软件 ， 进 行 样品 组 间 的 差异 表达 分 析 ， 将 差异 倍数 
llog2FoldChange|*1 H. padj«0.05 作为 筛选 差异 表达 基因 (DEGs) 的 条 件 。 对 符合 条 件 的 差 
异 表达 基因 进行 GO Ontology 和 KEGG Pathway 分 类 ， 利 用 超 几 何 检 验 对 注释 到 的 不 同 的 
GOCP<0.05) 和 KEGG 条 目 (P < 0.05 ) 进 行 富 集 分 析 。 采 用 FPKM 值 (Fragments per kilobase 
of transcriptper million fragments mapped) 来 展现 基因 表达 水 平 的 高 低 ， 在 R 软件 中 绘制 差 
异 基因 表达 量 的 火山 图 。 
2 结果 与 分 析 
2.1 转录 组 测序 与 组 装 

为 进一步 理解 不 同 品种 火龙 果 的 分 子 差异 机 制 ， 对 2 个 品种 火龙 果 共 6 个 样本 进行 

RNA-seq， 获 得 平均 7.71 Gb clean reads, Q30 均 在 90% 以 上 ， 总 碱 基 中 的 GC 数量 比例 为 
49.06%〈 表 1) ， 表 明 测 序 结果 良好 ，clean reads 质 量 合格 ， 可 进行 下 一 步 数据 分 析 。 
表 1 测序 数据 质量 表 
Table 1 Quality table of sequencing data 


样品 原始 Reads ”高 质量 Reads 数据 量 准确 率 >0.999 GC 含量 
Sample Raw Reads Clean Reads Clean bases Q30 (96) GC content (96) 
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BR_1 53 250 692 52 512 370 7.88Gb 91.36 48.98 
BR 2 50 237 344 49 478 466 7.42Gb 90.77 49.56 
BR 3 50 441 754 49 424 658 7.41Gb 90.14 48.69 
EY 1 58 014 828 56 782 758 8.52Gb 91.12 48.82 
EY 2 53 025 282 51 813 250 7.77Gb 90.99 49.16 
EY 3 49 516 914 48 245 168 7.24Gb 90.75 49.19 


ik: BRI, BR 2. BR 3 分 别 为 火龙 果 “ 白 肉 ” 的 3 个 重复 ，EY 1. EY 2. EY 3 分 别 为 火龙 果 “ 黄 龙 ” 的 3 个 重复 。 
Notes: BR. 1, BR 2 and BR 3 are three biological duplicates of pitaya ‘Bairou’; EY 1, EY 2 and EY 3 are three biological 


duplicates of pitaya ‘Ecuador Yellow’. 
2.2 火龙 果 转 录 组 Unigene 的 功能 注释 
使 用 BLAST 将 所 有 Unigene 5 NR, GO, KO, SwissPort, PFAM, KOG, NT 共 7 个 数 
据 库 进行 一 致 性 比 对 分 析 ， 对 各 数据 库 注 释 的 Unigene 数目 进行 统计 ， 进而 获得 火龙 果 转 录 
组 Unigene 的 功能 注释 信息 。 结 果 表 明 ，28 551 条 Unigene (45.02%) 在 NR 数据 库 比 对 成 
功 得 到 注释 ,在 GO, KO, SwissPort 等 数据 库 获 得 注释 的 Unigene 数目 依次 为 23 082(36.4%)、 
11 011 (17.3696), 21352 (33.6794) (图 2)。4 493 条 Unigene 同时 在 所 有 数据 库 中 注释 ， 
至 少 有 1 种 数据 库 注释 成 功 共 34 692 条 (54.71%), 28 714 条 未 获得 注释 。 
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Annotation of seven databases 
A. 注释 到 NR 数据 库 ，B. 注释 到 NT 数据 库 ，C. 注释 到 KO 数据 库 ，D. 注释 到 SwissProt 数据 库 ，E. 注释 到 PFAM 
数据 库 ，F. 注释 到 GO 数据 库 ，G. 注释 到 KOG 数据 库 。 
A. Annotated in NR; B. Annotated in NT; C. Annotated in KO; D. Annotated in SwissProt; E. Annotated in PFAM; F. Annotated 


in GO; G. Annotated in KOG. 
图 2 注释 结果 柱状 图 
Fig. 2 Histogram of Unigene annotation results 

2.3 不 同 火龙 果品 种 差异 表达 基因 分 析 

对 “白肉 ”和 “黄龙 ”火龙 果 的 差异 表达 基因 进行 分 析 。 结 果 表明 ，“ 白 肉 ” 和 “黄龙 
火龙 果 之 间 有 14 248 个 差异 基因 (padj < 0.05 H[log2FoldChange-10, B, “白肉 ”相对 
于 “黄龙 "， 有 5 446 个 基因 上 调 ， 表 达 量 高 于 “黄龙 "， 有 8 802 个 基因 下 调 ， 表 达 量 低 于 

“黄龙 ”( 图 3)。 
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© X, Y 坐标 轴 均 取 筛 选 条 件 值 的 对 数 形 式 ， 图 中 每 个 点 代表 一 个 差异 表达 基因 ， 红 色 圆 点 表示 有 显著 性 差异 的 上 调 基因 ， 
CN 绿色 圆 点 表示 有 显著 性 差异 的 下 调 基因 ， 蓝 色 圆 点 表示 无 显著 性 差异 的 基因 。FoldChange 表示 差异 倍数 ，padj 表示 校正 
E 后 的 pvalue。 
Y 
N X-axis and Y-axis present threshold value in log transform. Each dot is a differential expressed genes. Dots in red mean 
ed up-regulation genes, Dots in green mean down-regulation gene and dots in blue means non-regulation genes. FoldChange 
CN 
mu represents fold change, padj represents after correction pvalue. 


图 3 BR vs EY 间 基 因 差 异 表达 分 析 火 山 图 
Fig.3 Volcano graph of Differentially expressed genes in BR vs EY 
2.4 差异 表达 基因 GO 功能 分 类 与 富 集 分 析 

为 了 分 析 “白肉 ”与 “黄龙 ”差异 基因 的 生物 学 功能 ， 将 14 248 个 差异 基因 注释 到 GO 
Ontology 数据 库 中 ,分 别 从 生物 学 过 程 (biological process, BP )、 细 胞 组 分 (cellular component, 
CC) 和 分 子 功能 Cmolecular function; MF) 3 个 方面 ， 其 中 在 生物 学 过 程 分 类 中 ， 参 与 氧 
化 还 原 过 程 、 碳 水 化 合 物 代 谢 过 程 的 差异 基因 数目 最 多 , 分 别 为 632 和 371 个 ; 在 细胞 组 分 
分 类 中 ， 主 要 是 被 膜 等 的 差异 基因 组 多 ， 有 234 个 ; 在 分 子 功能 分 类 中 ， 主 要 是 氧化 还 原 酶 
最 多 ， 有 618 个 。 

UER PEDF 0.05 作为 富 集 标 准 。 在 生物 过 程 分 类 中 富 集 了 包括 DNA 整合 mRNA 
修饰 、mRNA 甲 基 化 、 有 丝 分 裂 、RNA 甲 基 化 等 39 个 条 目 ; 在 细胞 组 分 分 类 中 显著 富 集 了 
包括 肌 质 、 肌 质 网 等 28 个 条 目 ; 在 分 子 功能 分 类 中 显著 富 集 了 包括 四 吡咯 结合 、 血 红 素 结 
合 、 铁 离子 结合 等 54 项 。 结 果 说 明 ,“ 白 肉 ” 与 “黄龙 ”在 生长 发 育 、 细 胞 代谢 过 程 等 方面 
存在 差异 。 对 上 调和 下 调 的 差异 基因 分 别 进行 功能 分 类 ， 在 “白肉 ”与 “黄龙 ”差异 表达 基 
因 中 满足 P 值 小 于 0.05 的 GO 条 目 即 为 GO 功能 显著 富 集 ， 部 分 结果 见 图 4。 对 5 446 个 上 
调 基 因 进 行 GO 功能 分 类 与 富 集 分 析 ， 结 果 显 示 ， 在 生物 学 过 程 分 类 中 ， 参 与 氧化 还 原 过 程 

Coxidation-reduction process) 的 GO term 显著 富 集 ， 其 中 包含 差异 基因 341 个 ; 在 分 子 功 
能 分 类 中 ， 参 与 氧化 还 原 酶 活性 (oxidoreductase activity) 的 GO term 显著 富 集 ， 其 中 包含 
差异 基因 349 个 ， 表 明 “ 白 肉 ” 与 “黄龙 ”在 生长 发 育 、 细 胞 代谢 过 程 等 方面 存在 差异 ， 在 


Me 


分 子 功能 分 类 中 ， 抗 氧化 活性 Cantioxidant activity? 的 GO term 富 集 ， 表 明 “ 白 肉 ” 和 “ 黄 
Jb' 在 抗旱 性 方面 有 差异 。 在 下 调 基因 中 ， 主 要 集中 在 DNA 整合 (DNA integration) 等 方 
面 ， 差 异 无 显著 表现 。 
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Biological process Cellular component Molecular function 
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A: A， 氧 化 还 原 过 程 ，B.， 氧 化 还 原 酶 活性 ;，C， 抗 氧化 活性 ; D， 氧 化 还 原 酶 活性 ， 作 用 于 过 氧化 物 作 为 受 体 ，E， 铁 离 
子 结合 ，F， 四 吡咯 结合 ，6， 血 红 素 结合 。 


A.oxidation-reduction process; B. oxidoreductase activitys; C. antioxidant activity; D. oxidoreductase activity, acting on 


peroxide as acceptor; E. iron ion binding; F. tetrapyrrole binding; G. heme binding. 

B: A. DNA 整合 ，B. mRNA 修饰 ，C. mRNA 甲 基 化 ; D. RNA 甲 基 化 ， E. 病毒 沉 体 ; F. virion 部 件 ，G. RNA 定向 RNA 
聚合 酶 复合 物 ，H. mRNA 甲 基 转移 酶 活性 ; L 氧 裂解 酶 活性 ， 作 用 于 磷酸 盐 , J. AAE; K. ADP 结合 ; L. RNA 
甲 基 转 移 酶 活性 ，M. 腺 昔 甲 硫 氮 酸 依赖 性 甲 基 转移 酶 活性 ，N.RNA 定向 RNA 聚合 酶 活性 ; 

A. DNA integration; B. mRNA modification; C. mRNA methylation; D. RNA methylation; E. viral capsid; F. virion part; G. 


~ 


RNA-directed RNA polymerase complex; H. mRNA methyltransferase activity; I. carbon-oxygen lyase activity, acting on 
phosphates; J.  terpene synthase activity, K. ADP binding L. RNA methyltransferase activity, M. 
S-adenosylmethionine-dependent methyltransferase activity; N. RNA-directed RNA polymerase activity. 
图 4BR VS EY 上 调 (A) 下 调 B) 差异 表达 基因 的 GO 分 类 及 富 集 分 析 
Fig. 4 GO classification and enrichment analysis of up(A)-down(B)- regulation differentially 


expressed genes in BR VS EY 

2.5 差异 表达 基因 KEGG Pathway 富 集 分 析 

“白肉 ”与 “黄龙 ”共有 2 237 个 差异 表达 基因 注释 到 特定 的 KEGG 通 路 中 ,这些 通 路 主 
要 为 新 陈 代谢 、 生 物 合 成 、 糖 酵 解 、 信 号 传导 、 和 蛋白 质 加 工 等 ， 其 中 注释 到 为 新 陈 代谢 的 差 
异 基因 838， 占 差异 基因 总 注释 数目 的 37.5%， 注 释 到 生物 合成 的 差异 基因 335 个 ， 占 差异 基 
因 总 注释 数目 的 15.0%。 以 P< 0.05 作 为 篇 选 显著 富 集 的 条 件 ， “白肉 ”中 表达 量 较 “黄龙 
上 调 的 差异 林 因 主要 富 集 在 生物 合成 方面 ,例如 角质 、 木 栓 质 和 蜡 质 生物 合成 (cutin, suberine 
and wax biosynthesis) 、 茶 丙 素 的 生物 合成 (phenylpropanoid biosynthesis) ~ HIREA 
生物 合成 (indole alkaloid biosynthesis) 等 (图 5) 。 其 中 参与 角质 、 木 栓 质 和 蜡 质 生物 合成 
的 差异 基因 有 12 个 ， 分 别 是 HHT1、CYP86B1、CYP86A1、FAR、CYP86A1、HHT1、PXG、 
CYP86A4S、K15404、CER1 FAR K15404, CERI K15404、CER1。“ 白 肉 ” 中 表达 量 较 “ 黄 
JE' 下 调 的 差异 基因 主要 富 集 在 生物 合成 方面 ， 例 如 类 倍 半 莫 和 三 熙 的 生物 合成 
sesquiterpenoid and triterpenoid biosynthesis、 精 氢 酸 的 生物 合成 Carginine biosynthesis) 、 脂 
肪 酸 的 生物 合成 (fatty acid biosynthesis) 等 (图 6) 。 


富 集 因子 Rich factor 


A. 维生素 B6 代 谢 ; B. 泛 柄 和 其 他 布 类 醒 生 物 合成 ; C. 栈 氨 酸 代谢 ; D. RNA RS: E 喀 喧 代谢, F EARE: G. 光 
合作 用 天 线 蛋 白 ; H 光合 作用 ; I 茶 丙 素 的 生物 合成 ，J. 茶 丙 氨 酸 代谢 ，K. 戊 糖 和 葡 糖 醛 酸 相互 转化 ，L. 氧化 磷酸 
化 ; M. 嘻 哄 生物 碱 生 物 合成 ，N. 甘氨酸 、 丝 氨 酸 和 苏 氨 酸 代谢 ，O. 半 胱 氨 酸 和 和 蛋氨酸 代谢 ;，P. 角质 、 木 栓 质 和 蜡 质 
EWER: Q. 甜菜 色素 生物 合成 ，R. 抗坏血酸 盐 和 醛 酸 盐 代谢 ，S. 精 氨 酸 和 膊 氨 酸 代谢 ，T. 氨基 糖 和 核 苷 酸 糖 代 谢 。 


A. Vitamin B6 metabolism; B. Ubiquinone and other terpenoid-quinone biosynthesis; C. Tyrosine metabolism; D. RNA 


ChinaXiv 合 作 期 刊 
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polymerase; E. Pyrimidine metabolism; F. Proteasome; G. Photosynthesis - antenna proteins; H. Photosynthesis; I. 
Phenylpropanoid biosynthesis; J. Phenylalanine metabolism; K. Pentose and glucuronate interconversions; L. Oxidative 
phosphorylation; M. Indole alkaloid biosynthesis; N. Glycine, serine and threonine metabolism; O. Cysteine and methionine 
metabolism; P. Cutin, suberine and wax biosynthesis; Q. Betalain biosynthesis; R. Ascorbate and aldarate metabolism; S. 
Arginine and proline metabolism; T. Amino sugar and nucleotide sugar metabolism. 
图 5 “白肉 ”与 “黄龙 ”火龙 果 差 异 表达 上 调 基因 富 集 前 20 的 KEGG 通 路 
Fig.5 BR VS EY Differentially expressed down-regulated genes were significantly enriched in 
the top 20 KEGG pathways 
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A. AR TARMA REAREA B. 类 倍 半身 和 三 熙 的 生物 合成 ; C. RNA 降 解 ，D. 核 黄 素 代谢 ，E. 丙酮 酸 代谢 ; 
F. 茶 丙 素 生 物 合成 ，G. 过 氧 物 酶 体 ; H. 其 他 类 型 的 O- 聚 糖 生物 合成 ; 工 其 他 聚 糖 降解 ， J. 糖 基 碰 脂 酰 肌 醇 - 锚 生物 合 
Wü. K. 糖 酵 解 / 糖 异 生 ，L. 果糖 和 甘露 糖 代谢 ，M. 脂肪 酸 降解 ，N. 脂肪 酸 的 生物 合成 ，O. BRIE- P. 类 胡 
萝卜 素 生物 合成 ，Q. 光合 生物 的 固 碳 作用 ; R 生物 素 代 谢 ，S. 精 氨 酸 的 生物 合成 ，T. a -亚麻 酸 代谢 。 


A. Valine, leucine and isoleucine degradation; B. Sesquiterpenoid and triterpenoid biosynthesis; C. RNA degradation; D. 


Riboflavin metabolism; E. Pyruvate metabolism; F. Phenylpropanoid biosynthesis; G. Peroxisome; H. Other types of O-glycan 
biosynthesis; I. Other glycan degradation; J. Glycosylphosphatidylinositol(GPI)-anchor biosynthesis; K. Glycolysis / 
Gluconeogenesis; L. Fructose and mannose metabolism; M. Fatty acid degradation; N. Fatty acid biosynthesis; O. Circadian 
rhythm - plant; P. Carotenoid biosynthesis; Q. Carbon fixation in photosynthetic organisms; R. Biotin metabolism; S. Arginine 
biosynthesis; T. alpha-Linolenic acid metabolism. 
图 6 “白肉 ”与 “黄龙 ”差异 表达 下 调 基因 富 集 前 20 的 KEGG 通 路 
Fig.6 BRVS EY Differentially expressed up-regulated genes were enriched in the top 20 KEGG 
pathways 


3 讨论 
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KE (张驰 等 ，2015) ， 桂 花 〈 张 雪松 等 ，2016) . HER CERA, 20190 等 。 本 研究 中 


我 们 利用 RNA-seq 技 术 对 火龙 果 不 同 品种 进行 转录 组 学 分 析 ， 基 于 良好 的 测序 结果 ， 我 们 进 
步 分 析 “ 白 肉 ” 与 “黄龙 ”火龙 果 二 者 之 间 的 差异 表达 基因 ， 发 现 共有 14 248 个 差异 表达 
基因 ， 通 过 GO Ontology 分 析 表 明 这 些 差异 基因 主要 富 集 在 新 陈 代谢 、 生 物 合成 、 糖 酵 解 、 
信号 传导 、 和 蛋白 质 加 工 等 过 程 ， 表 明 “ 黄 龙 ” 可 能 在 生长 发 育 、 抗 性 方面 与 “白肉 ”有 差异 。 
在 分 子 功能 分 类 中 ， 相 较 于 “黄龙 ”， “白肉 ”火龙 果 上 调 表 达 的 基因 中 ， 参 与 氧化 还 原 酷 
活性 Coxidoreductase activity) 的 差异 基因 349 人 个， 表明“ 白肉 ”与 “黄龙 ”在 生长 发 育 、 细 
胞 代谢 过 程 差异 显著 ，“ 白 肉 ” 比 “黄龙 ”的 生长 发 育 更 快 ， 代 谢 更 加 旺盛 ， 与 实际 栽培 
植 过 程 中 的 现象 相符 ，“ 和 白肉” 表现 出 比 “ 黄 龙 ” 更 快 的 生长 速度 ， 在 一 年 中 结果 的 批 次 ， 
果实 的 发 育 都 有 差异 ，““ 和 白肉” 一 年 结果 大 概 8 批 ， 从 授粉 到 成 熟 需要 28~30 d，“ 黄 龙 ” 一 
年 结果 2 批 , 夏 果 和 冬 果 , 成 熟 期 需要 90~120 d. 在 分 子 功 能 分 类 中 , 抗 氧化 活性 Cantioxidant 
activity) 差异 显著 , 表明 “白肉 ”和 “黄龙 ”在 抗旱 性 方面 有 差异 。 为 了 在 干旱 环境 中 生存 ， 
仙人 掌 科 植物 必须 具备 干旱 反应 机 制 , 特别 是 抗 氧化 防御 系统 。 持续 的 干旱 胁迫 会 导致 氧化 
性 物质 如 QO?、H?0*、0”-、OH 的 积累 ， 这 些 物 质 会 损害 细胞 ， 甚 至 导致 死亡 (Fang & Xiong, 
20150 ， 这 些 抗 氧 化 活性 基因 在 “白肉 ”中 富 集 ， 上 调 表 达 ， 可 能 使 得 “白肉 ” 抗 干旱 能 

强 于 “黄龙 ”。 
在 KEGG Pathway 中 以 已 < 0.05 作为 富 集 条 件 ， 在 “白肉 ”上 调 表达 的 基因 中 ， 主 要 富 
集 在 角质 、 木 栓 质 和 晴 质 生物 合成 ， 其 中 参与 角质 、 木 栓 质 和 量 质 生物 合成 的 差异 基因 有 
12 个 ， 主 要 代表 有 CYP86 和 CER1。 植 物 蜡 质 是 植物 在 进化 过 程 中 为 了 适应 环境 的 变化 而 
演变 出 来 的 一 种 结构 ， 表 皮 旺 质 在 植物 抵御 外 界 不 良 环境 、 保 持 组 织 和 器 官 功 能 、 保 证 植物 
正常 发 育 等 方面 起 了 重要 的 作用 〈( 胡 晓 敏 等 ，2007;， 李 魏 强 等 ，2006; 任 春 涛 等 ，2019) 。 
上 表皮 蜡 沉 积 的 主要 功能 是 减少 表皮 的 水 分 流失 ,这 一 特性 有 助 于 增强 耐 旱 能 力 。 此 外 , 在 
与 草食 性 昆虫 和 植物 病原 真菌 的 相互 作用 中 起 着 重要 作用 (Aarts et al., 1995) , WE EE 
定 程 度 上 可 以 使 植物 免 受 低温 胁迫 作用 〈 李 婧 婧 等 ，2011) 。 其 中 12 个 关键 基因 参与 角质 、 
亚 黄 素 和 蜡 的 生物 合成 ， 如 CYP86 编码 @- 基 脂肪 酸 的 形成 ， 促 进 在 软木 脂 生物 合成 ， 形 成 
角质 层 。 这 些 脂肪 类 生物 聚合 物 可 以 防止 植物 水 分 和 营养 流失 以 及 病原 体 的 感染 (Aarts et al., 
1995; Hofer etal., 2008) . CERI 编码 参与 长 链 醛 转化 为 烷烃 的 蛋白 ， 这 是 植物 体 中 螨 生 
物 合成 的 关键 步骤 (Hofer etal., 20080 。 本 研究 对 这 些 关 键 调控 基因 分 析 发 现 ， 此 类 基因 在 
“白肉 ”中 具有 较 高 的 表达 量 ， 这 些 可 能 与 植物 的 抗旱 和 抗 病 有 关 。 
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